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j^jfj, Abstract 

Chaotic transients occur in many experiments including those in fluids, in simulations of the plane 
Couette flow, and in coupled map lattices and they are a common phenomena in dynamical systems. 
Superlong chaotic transients are caused by the presence of chaotic saddles whose stable sets have fractal 
dimensions that are close to phase-space dimension. For many physical systems chaotic saddles have a 
big impact on laboratory measurements, and it is important to compute the dimension of such stable 
| sets including fractal basin boundaries through a direct method. In this work, we present a new method 

(J , to compute the dimension of stable sets of chaotic saddles directly, fast, and reliable. 

Chaotic transients are common phenomena in dynamical systems. They were observed in the Lorenz 
system ^ |2j, in fluid experiments El Eli m many low-dimensional dynamical systems El HI El- i n 
spatiotemporal chaotic dynamical systems ^3 E]i simulations of the plane Couette flow (a shear flow 
between two parallel walls) E21E|, and in coupled map lattices ^3- Generally, these transients are caused 
by the presence of chaotic saddles ^3 HH1 El EE] (i-e., strange repellers). When there is a chaotic saddle 
in the phase space, trajectories originating from random initial conditions in a neighborhood of the stable 
set of the chaotic saddle (such as in the neighborhood of a fractal basin boundary (BB)), usually wander 
in the vicinity of the chaotic saddle for a finite amount of time before escaping the neighborhood of the 
chaotic saddle and settling into a final attractor. Chaotic saddles can strongly effect observed experiments 
and a "turbulent" state in experiments could be supported by a chaotic saddle rather than an attractor. 
\& • If there are two or more attractors in the system, then chaotic saddles may be buried in the boundaries 

of basins of attraction leading to fractal BBs. Supertransients are superlong (chaotic) transients and occur 
commonly in spatiotemporal chaotic dynamical systems ^JJEIE]- I n such a case, trajectories starting 
from random initial conditions in a neighborhood of the stable set of a chaotic saddle wander chaotically for 
a long time before settling into a final attractor. It was observed in numerical experiments that spatially 
extended systems exhibit supertransients so that the observation of the system's asymptotic attractor is 
practically impossible ^JEJ- Reference ^I] investigates the geometric properties of the chaotic saddles 
which are responsible for the supertransients in spatiotemporal chaotic systems, and reports that the stable 
set of the chaotic saddle possesses a fractal dimension which is close to the phase-space dimension. There 
are a few tools for computing the dimension of (one-dimensionally unstable) stable sets of chaotic saddles 
[§1 1191 ITK1 1281 1171 0(11 II lj . These methods involve the computation of Lyapunov exponents. A frequently used 
method is the combination of the Sprinkle method !) and the method for computing the maximal Lyapunov 
exponent of a typical Saddle Straddle Trajectory ^3I231E|- However, computing the Lyapunov exponents 
of a Saddle Straddle Trajectory in a supertransient system is extremely time consuming. 

The purpose of this Letter is to present a simple, general method (involving no Lyapunov exponents) 
for computing the dimension of the stable set of a chaotic saddle that is responsible for the supertransients 
and the computation is direct, fast and reliable. We refer to this algorithm as the "Straddle Grid Dimension 
Algorithm" (SGDA); it can be applied to flows (differential equations) or maps. Our main result is that 
the SGDA can be applied to supertransient systems having one-dimensionally unstable stable sets such as 
basin boundaries in d-dimensional space (d > 1). To illustrate the effectiveness of our method for computing 
dimension see Fig. El As indicated by our numerics, the SGDA is superior to applying the box-counting 
algorithm (BCA) by utilizing straddle pairs for the end points of the boxes. 
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We present the SGDA for maps but only trivial modifications are necessary for flows 1 . Let / be a 
continuous map from the d-dimensional phase space M into itself such that / is a C 3 -diffeomorphism when 
d > 2, or a C 2 -map when d = 1. Let R be a "transient region" in M (that is, a region in the phase space that 
contains no attractor; of course, a transient region may contain points that belong to an attractor). Under 
some mild assumptions, almost all orbits with initial conditions in the transient region leave this region in a 
finite number of iterates under / and approach one of the attractors of /. The escape time T(x) of x from 
R is the minimum integer n > for which the nth iterate of x is not in i?; write T(x) — oo if all forward 
iterates of x are in R. The escape time of points in the transient region is at least 1, so for x G R, one has 
T(x) > 1. Let C be the largest invariant set of / in R, that is, /(C) = C, so T(x) = oo for x G C. Assume 
that C is nonempty and T(x) is finite for almost every point x G R. The maximal invariant C can contain 
complicated chaotic sets but it can be just a periodic orbit or a fixed point. The stable set of C, S c , is the 
collection of all points x G R such that T(x) — oo. If the escape time T(x) of a point x is high, then the 
initial condition x is close to S c . 

Before we present our method, we first discuss the usual indirect method for computing the dimension 
of stable sets. Assume that C is a chaotic saddle in the transient region R. To see what happens for a 
uniform distribution of initial points, assume that at t = 0, one selects N(0) initial points (for N(0) large) 
from a uniform distribution on R. Evolve these iV(0) initial conditions under the dynamics. For t > 0, 
let N(t) denote the number of trajectories that are still in the region R at time t. Typically, N(t) decays 
exponentially as time t increases |2T]. Let A t be the set of points x G R with escape time T(x) > t. Write 
m(At) for the Lebesgue measure of A t . Typically, m(A t ) — » exponentially fast as t — » oo; that is, for 
large t, m(A t ) K exp(— t/r), for some positive constants K and t; (t is called the average escape time 2 ). 
For two-dimensional systems , the dimension d u of the collection of points of intersection of a line segment 
with the stable set S c can be computed indirectly from the relation r « 1/(1 — d u )Xi OEIj where Ai is 
the maximum Lyapunov exponent for typical trajectories on the chaotic saddle C. By using the Sprinkle 
method 3 , one computes r, and by using the PIM-triple method 4 , one computes Ai. Thus, it is important 
to have a general method that computes such dimensions directly, fast and reliable. 

We now introduce some concepts needed for describing our new method. For points a, b G M, and e > 0, 
we call {a, 6} an e-pair of points if the distance \\a — b\\ = e, and a straddle pair if the line segment [a, b] 
contains points of S c . From the computational point of view, to decide if a pair {a, 6} is a straddle pair, 
one performs a finite (and usually small) number of computations studying a refinement of the pair {a, b}. 
A S-refinement of the pair {a, b} is a finite set of points X = {Xi}f =0 in the segment [a, b] with Xq = a 
and X s = b and the distance between two consequtive points Xi and Xi + \ is S\\a — b\\. In our method for 
computing the dimension of stable sets of a chaotic saddle, we make use of a test (based on the PIM-triple 
concept ^H1E3E|)> for verifying whether a pair {a, b} is a straddle pair. A triple {Xi,Xj,Xk} of points 
in a typical line segment in M intersecting S c is called a PIM triple if (a) Xj is between Xi and Xk, and 
(b) the escape time of Xj is bigger than both the escape times of Xj and Xk, that is, T(Xj) > T(Xi) and 
T(Xj) > T{X k ). The "Straddle test" below is a basic tool for the SGDA. 

Straddle test. Let / denote a typical line segment in M intersecting S c and e > be given. Let {a, 6} 
be an e-pair in / and X — {Xi}f =Q be a refinement of the pair {a, 6} in /. If X = {Xi}f =0 includes a PIM 
triple then {a, b} is a straddle pair 5 . 

Note that the above procedure tests only finite sets of points; even if the Straddle test does computa- 

x To apply the method to flows, consider the map generated by the numerical solver which takes the system from a point in 

the phase space at time t n to a point at time tn+l = t n + At„, and apply the method. 

2 An average escape time t of the chaotic transients can be defined by X/r = lim t ^oc lim^o)^ ± In GJH3E]- 
3 The Sprinkle method starts with N initial points x that are uniformly distributed in the transient region. Given two integers 

Ni and Nf, for each initial point x with T(x) > (JVj + AT*), the first JVj iterates of x are discarded and so are the last Nf, 

the remaining iterates are plotted. The accuracy and the speed of the Sprinkle method for finding points close to the chaotic 

saddle are limited 0. 

4 Unlike the Sprinkle method, the PIM-triple method quickly generates trajectories which are very close to the stable set of 
the chaotic saddle, and Lyapunov exponents can be computed fast and accurately for non-supertransient systems 1 1 "il I2.j| 1 1 71 . 

5 If the stable set in R is a BB (that is, a point x £ M is on the BB if every neighborhood of x contains points of at least 
two different basins), then instead of the above straddle test, we may choose to apply the following "BB straddle test". BB 
straddle test. Let / denote a typical line segment in M intersecting S c and e > be given. Let {a, b} be an e-pair in J and 
X = {Xi}f_ Q be a refinement of the pair {a, b} in I. If there exist < i,j < s such that Xi and Xj belong to different basins, 
then the pair {a, b} is a straddle pair. 
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Figure 1: The left refinements (grids) used in the selection procedure. The points of the refinements are 
marked by short vertical segments. 
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Figure 2: (a) \D g — D\ versus D (black) and \Db — D\ versus D (gray), where D g (Db) is the numericaly 
calculated grid (box) dimension, and D is the dimension calculated analytically for our first numerical 
example. We used n = 5, k = 17, and r — 3; (b) D g versus [i (black) and Db versus /j, (gray), where D g (Db) 
is the numericaly calculated grid (box) dimension, and \x is the parameter of our second numerical example. 
We used n = 5, k = 17, and r — 3. The refinements in (a, b) are tested with the Straddle test. Our numerics 
clearly illustrate a precision gain up to an order of magnitude when the dimension is calculated by our 
method as compared to when the dimension is calculated by the box counting algorithm. Both calculations 
use the same grid of points so they are equally fast. Indirect methods involving the calculation of Lyapunov 
exponents are much slower; (c) Similar to (b) for the last numerical example except that n = 5, k = 19, and 
r = 1. The refinements are tested with the BB straddle test. 



tionally diagnose a pair {a, b} as not a straddle pair for a given refinement, this does not imply that the pair 
{a, b} is not a straddle pair. Quite possibly, a different choice of the refinement of {a, b} may reveal that 
{a, b} is a straddle pair. Thus, computationally, one always underestimates the number of straddle pairs. 
The BCA has specific procedures for selecting line intervals of size e in /. Then, the ends of the intervals 
are chosen as e-pairs, and each e-pair may be further refined. The BCA evaluates the fraction of e-pairs 
that are straddle pairs according to the Straddle test, where e decreases on a logarithmic scale. The fraction 
of straddle e-pairs detected computationally is used as an estimate of the fraction of line intervals of size e 
that straddle the fractal set. Thus, the fact that every test based on a finite number of computations does 
not accurately detect straddle pairs is a major computational difficulty for the BCA. We now present the 
SGDA which allows us to compute the dimension of the stable set of a chaotic saddle. Given a grid of points 
in /, we use the grid points to construct (5-refined e-pairs that we test with the Straddle test. We denote 
the fraction of the selected pairs that are straddle pairs by ip e . Similar to the scaling of fractal and uncer- 
tainty dimensions |22l I16| , using an approach similar to |lf>| guaranteeing that these types of dimensions of 
certain basin boundaries coincide, we argue that ip e ~ e 1_ '° 9 for small e 6 , where D g is by definition the 
straddle grid dimension. An important ingredient for our algorithm is the selection procedure for e-pairs 
and ^-refinements. In the "Selection procedure" below, we consider left ^-refinements of {a, 6}; that is, a 

6 The probability that an interval of size e contains a straddle pair of any size less than e must equal the probability 
that the interval straddles the fractal set. Using Lcbcsgue measure to count the number of (5-refinements, we have </> € = 
[J"q (e — S)ijigdS\ I [Jg(e — S)d5j = (2/e 2 ) J S (€ — 5)ip$d5. Taking two derivatives of e 2 rj> e with respect to e we obtain ip e = 
<j>c + ed e <f>e + e 2 fe0e/2, and thus, since e ~ e 1 " 15 so does ip e . 
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(5-refinement of {a, b} in which the right point b is excluded. Let a and b be the end points of /. For fixed 
nonnegative integers n, k, and r, start with a grid G n + k + r which is a left 2-("+ fc+r )-refinement of {a, b}. 
Before we state the selection procedure, we specify subgrids to be used. G n is the left 2 _ ™-refinement of 
{a, b} included in G n + k + r . Denote the grid points of G n by where Y\ = a and Yi < Y. L+ i for all 

i. For all i, write Gfc+ r [i] for the left 2 ~( fc +r) -refinement of {Yi,Yi + i}, G r [i] for the left 2 _r -refinement of 
{Yi,Y i+ i} (see Fig. EJ, and S^(G r [i]) for the resulting subgrid of Gk+ r [i] when G r [i] is shifted j grid steps 
to the right. The union of the S j {G r [i\Ys (0 < j < 2 k ) is G k+r [i]. 
Selection procedure for e-pairs and 5 -refinements. 

1. Start with the left 2~"-rcfinemcnt of {a, b} which is the grid G„. Initialize i = 1. (Of course, here 
e = 2-™.) 

2. Start with the 2~ n -pair {Yi,Y i+ i}, the left 2~ r -refinement of {Yi,Y i+ i} (grid G r [i]), and the left 
2-( fc + r )-refmement of {Y % , Y i+1 } (grid G k+r [i\); see Fig.HJ 

3. Select the subgrids S j (G r [i]) in the grid G k+r [i] for < j < 2 k . 

4. Increase i by one. Repeat steps 2., 3., and 4. for all 1 < i < 2™. 

When shifting the grid, the resulting grid S^(G r [i]) is not a refinement of {Yi, i^+i} but it is a left 2~ r - 
refinement of {S^(Yi), (Yi + i)}. The above procedure selects a total of 2 n+k refinements of 2 r points each. 
Every point of the grid G n + k +r is considered exactly once in the selected refinements. 

Straddle Grid Dimension Algorithm. 7 

(i) For three integers n, k, and r, consider a grid G fl+k - +r which is a left 2-( ft +^+'')-rcfinemcnt of {a, b] 
where a and b are the end points of /. Initialize k = k and n = h\ 

(ii) For each point x € G n ,^, r compute the escape time T(x) from the region R; 

(iii) Apply the selection procedure for n, k, and r. Using the Straddle test, count the number of the 
selected refinements whose pairs straddle S c . Denote the ratio of this number of selected refinements to the 
total number of refinements by ip n (= ip e where e = 2~"); 

(iv) Increase n by 1, decrease k by 1. If k < proceed to (v), else proceed to (iii); 

(v) Plot log(-0„) versus n log 2. The slope of this graph is 1 — D g , where D g is the straddle grid dimension 
on /. 

Denote the fraction of line intervals of size 2~ n that straddle the fractal set S c by (j> n . The number of 
refinements tested by the SGDA at step (iii) is constant 2 n+k = 2 n+k , independent of k or n. The SGDA 
integrates ip n by a rectangle method using grid points, just as the BCA integrates (j> n . However, the SGDA 
uses all of the grid points at every step of the algorithm, only sorted differently, in order to get the i\> n 
estimates. The SGDA performs a total of k2 n+k refinement tests, to be compared with the number of tests 
of the BCA which is only 2™(2 fc ~ 1 - 1). Thus, the V™ values are significantly more accurate than the <f> n 
values. We note that with increasing n while keeping the same grid, the BCA performs better since values 
of 4> n with low n and high error are removed from the least square fit toward the fractal dimension. Under 
the same conditions, the SGDA is expected to perform worse. In this case, all ip n values have comparable 
error, and eliminating some of them from the least square fit increases the error of the fitted slope. 

Numerics, (a) A piecewise linear dynamical system. We consider the one-dimensional map |24| 

r {x + !)/ ii- l, x < -1/(1 - H)/{v + m); 
9v A x ) = \ -x/u, \x\ < v{\ - (i)/(v + n); (1) 

[ (x-l)/n + l, v{l - m)/(^ + M) < x; 

where < v < 1, and < /i < 1/2. For (v + 2/i) < 1, the dynamical system has two attractors at 
±00 that share a fractal BB. Here we choose v = /i. The BB points are included in [-1,1], and the BB has 
a unique dimension which is D = — log 3/ log /i. Fig. Efa) shows the box counting and the straddle grid 
dimension for 0.935 < D < 1 (i.e., 0.31 < /1 < 0.3325) where we used the Straddle test for refinements. The 
black line represents \D g — D\ versus D, and the gray line represents \Db — D\ versus D, where is the 

7 A fractal BB is a one-dimensionally unstable stable set of a chaotic saddle. For computing the dimension of fractal BB, 
instead of the box-counting method, one may apply the Uncertainty Dimension Algorithm (UDA) 1221 . If the stable set is a 
BB, one may apply the BB straddle test instead of the Straddle test. If / intersects a BB, then step (ii) in the SGDA may be 
replaced by "For each point x £ G n +fe+ r except the right point of / determine the basin the point x is in." and the first part of 
step (iii) may be replaced by "Apply the selection procedure for the integer values n, k, and r and then the BB straddle test." 
When the stable set is a fractal BB, SGDA is superior to UDA. 
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Figure 3: D g versus n (black) and Db versus /x (gray), where D g (Db) is the numericaly calculated grid (box) 
dimension, and /i is the forcing amplitude of the Duffing oscillator. We used n = 1, k = 18, and r = 1. 
The refinements are tested with the BB straddle test. Again, our numerics clearly illustrate a considerable 
precision gain of SGDA as compared to BCA. Both algorithms use the same grid of points so they are equally 
fast. A jump in dimension is expected at /i sa —0.2446 due to a saddle-node bifurcation. The graph of D g 
clearly illustrates this jump, while the graph of Db does not. 



box counting dimension. We see that the SGDA performs significantly better than the BCA using the same 
grid of points, (b) Nonlinear examples. We consider the logistic map f^{x) = /ix(l — x) for < x < 1 and 
fi > 4; it is well-known that this system has a chaotic saddle. Figure |2Jb) shows D g (black) and Db (gray) 
versus fi. The variation of D g has much less fluctuations than that of Db, suggesting that the computation 
of D g is more accurate. 

Our next example is constructed as follows. Start with f^(x) for fi = 3.832 where there is a stable period 
three orbit; denote its third iterate by F(x). Consider now the map /i M (x) = F(x) + fixe~ 7x . The map 
has fixed points attractors that share a fractal BB for < fi < 0.02. Figure EJc) shows D g (black) and Db 
(gray) versus \x. The accurate dimension graph generated by the SGDA method provides interesting and 
useful information on certain bifurcations of h^. For example, the jump in the dimension graph at /i « 0.01 
in Fig. EIc) is caused by a saddle node bifurcation. 

The last example that we present is obtained as follows. Consider the forced Duffing oscillator |25H2Kll77) 
given by d 2 x/dt 2 + 0.15 dx/dt + x — x 3 = /icos(i). Our dynamical system is the two-dimensional Poincare 
map at constant phase of the forcing t mod 2ir — ir. The map has attracting periodic orbits whose basins 
share a fractal boundary for —0.247 < /i < —0.242. In order to calculate the fractal dimension using SGDA, 
we choose a line segment / in the phase space of the variables x and dx/dt. All the points that belong to 
this segment have dx/dt = 1.35 with x ranging from —2 to 2. The numerical results obtained for the fractal 
dimension by applying both BCA and SGDA arc presented in Fig. [31 Again, our numerics clearly illustrate 
a considerable precision gain of SGDA as compared to BCA. Both algorithms use the same grid of points 
so they are equally fast. A jump in dimension is expected at fj, « —0.2446 due to a saddle-node bifurcation 
which is similar to the one mentioned in our previous example. The graph of D g (black) in Fig. 01 clearly 
illustrates this jump, while the graph of Db (gray) does not. 

Summary. We introduced a novel method for computing the dimension of both fractal BB and the 
stable sets of chaotic saddles in supertransient systems. The SGDA is the first general, but simple method 
for computing the dimension of stable sets of chaotic saddles directly, fast and reliable. We expect the SGDA 
to be a useful additional tool in analyzing the dynamics of fluids, lasers and transition to turbulence in shear 
flows. 
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